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I. INTRODUCTION 



The problem of quantum chaos in intrinsically degenerate systems possesses a number of 
interesting properties. The KAM theorem in these systems is not applicable^ and in certain 
models an arbitrary small perturbation is sufficient to induce an infinite stochastic web in 
phase space. The character of the web is determined by the type of perturbation. When 
the system is perturbed by a monochromatic wave, the web width exponentially decreases 
with increasing action / and the motion is practically localized.! If the perturbation has the 
form of periodic 5 - impulses, the web width is constant and the particle, traveling along 
the web, can diffuse to infinity! (see also Refs.i'@). In both cases the web in (p,x) phase 
space has a crystalline or quasicrystalline structure. The chaotic motion in these systems 
has been termed weak chaos since chaos occupies only a small portion of phase space ( See 
e.g., Ref. 2). In contrast, in strong ( or global ) chaos, the web structure disappears and 
most of the phase space is filled with chaotic orbits. The peculiar resonance structure and 
the appearance of an infinite stochastic web make degenerate systems very attractive objects 
in which to study quantum manifestations of chaos. 

The problem of weak and strong chaos has been mainly explored within the context of 
the kicked harmonic oscillator!!! It has been shown that the time of classical description 
of quantum averages is considerably longer for weak than for strong chaos. The role of the 
symmetry of the quasienergy (QE) functions was also analyzed in Ref.i. It was found that 
under certain conditions, quantum diffusion within the stochastic web was truncated by 
quantum interference effects!, similar to the case of strong chaos in the kicked rotorS. The 
problem of quantum chaos on the stochastic web has also been studied intensively in recent 
years within the framework of the generalized kicked Harper modelBi Another simple and 
yet important (especially to solid state physics) example of a degenerate system is that of a 
charged particle moving in a uniform magnetic field and interacting with a monochromatic 
wave, propagating perpendicularly to the magnetic field direction under the condition of 
cyclotron resonance. This 2D problem is tantamount to the ID harmonic oscillator in a wave 
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field. BI12I In our previous worked we studied this system quantum mechanically, focusing on 
the resonance approximation, which appears as the first order perturbation for the Floquet 
Hamiltonian. In the quasiclassical limit the dynamics in the resonance approximation is 
globally regular in phase space. The structure of the Floquet spectrum and the QE (quasi- 
energy) eigenf unctions for the exact and near resonance cases were obtained and related 
to the classical phase space structure. The evolution of various representative initial states 
was investigated and the close connection between classical and quantum dynamics at the 
cyclotron resonance was demonstrated. In Ref.@ it was shown that the boundaries of the 
quantum cells act as dynamical barriers to the probability flow. In the quasiclassical limit 
the dynamical barriers were found to correspond to the separatrices in classical phase space 
and tunneling through the "quantum separatrices " was explored numerically. 

In the present work we study the dynamical effects of chaos on the above mentioned 
system under the condition of cyclotron resonance u = uo c (u and u c are the wave and 
cyclotron frequencies, respectively). This model seems to be more closely related to exper- 
imental realizations in solid state physics than the kicked system. The evolution operator 
for an arbitrary number of periods of the external field is built and its structure is explored 
in terms of the QE eigenstates under the conditions of weak chaos. The structure of the 
evolution operator matrix is more complex than the more typical band-like matrix structure. 
Thus, the usual diagnostics of quantum chaos predicted by band random matrix theoryEiiE 
do not work in the case of weak chaos because only a small number of QE eigenstates is 
affected by the perturbation. A new phenomenon of diffusion over the quantum separatrices 
and the effect of weak chaos on the diffusion are investigated here and compared in the 
quasiclassical limit with the dynamics in classical phase space. The crossover from purely 
quantum diffusion to a diffusion which in the quasiclassical limit corresponds to classical 
diffusion within the stochastic web, is determined. The width of the quantum separatrices 
is computed and compared with the width of the classical stochastic web. 

The paper is organized as follows. In Sec. II the basic model is introduced. Also in Sec. 
II, the structures of the nonstationary Schrodinger equation and the evolution operator in Hq 
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representation are discussed. In Sec. Ill the properties of the QE eigenstates are described. 
The phenomenon of diffusion via the quantum separatrices and the influence of chaos on the 
diffusion are investigated numerically in Sec. IV. In Sec.V we draw our conclusions. 



II. THE EVOLUTION OPERATOR 

The Hamiltonian of a charged particle in a magnetic field interacting with a monochro- 
matic wave reads 

E= KV ^ ' +v cos(kx -cut) = H + V{x,t), (I) 

where m and e are, respectively, the mass and charge of the particle, p is the momentum, 
k is the wave vector, u is the wave frequency, and Vq is the amplitude of the perturbation. 
We choose the gauge of A in the form A = (0,Hr,0) so as to have the magnetic field 
H along the z-direction and to have the momentum p y as an integral of motion. The 
Hamiltonian is equivalent to a ID simple harmonic oscillator perturbed by a monochromatic 
wave field. Hence, the problem is to determine the dependence of the wave function on only 
two variables, x and t. 

It is convenient to expand the state vector in the harmonic oscillator basis 

if)(x, t) = J2 Cn{t)^ n {x) ex P (-iE n t/h), (2) 

n 

where ip n {x) is the nth eigenfunction of the simple harmonic oscillator Hamiltonian H and 
E n = hui c (n + 1/2) is the energy of the nth Landau level. Using (fj), the nonstationary 
Schrodinger equation 

ih^^=mx,t) (3) 

yields a set of differential- difference equations for the coefficients C n (t), 

ihC n = v J2 [V^l+m smM) + Kg +m cosM] C n+m e- m ^. (4) 
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The matrix elements V^ +m (V^„ +m ) describe the transitions between the levels of opposite 
(equal) parity and can be expressed via the Laguerre polinomials as0 



yd) j-l) m h m e 4 r 2m+lA , , 

v n,n+2m+l ~ ~ / \ 0d ) 

2 m+1 ^(n + l) . . .(n + 2m+l) 2 



Vn,n+2m — - r = = -LrT( n )i (5b) 

2 m+1 ^/(n + l) ...(n + 2m) 2 
where h = (ka) 2 plays the role of an effective (dimensionless) Planck constant and a = 



m 



hc/eH is the magnetic length. For n ^> 1 ^> h the matrix elements can be approximated 
terms of the Bessel functions J m of order m byi 

(1) i (-l) m n m+1/2 e~4 

K,n+2m+i = o T, \ = J 2m+ i{v2nh), (6a) 

2 J(n + 1) . . . (n + 2m + 1) 



vZ +2m = - , l ~ " e 4 J^)- (6b) 
2 J(n + 1) . . . (n + 2m) 

Since the perturbation is periodic in time, Floquet theory can be used to describe the 
time evolution of the system in terms of the QE spectra e q and the QE eigenfunctions 
ip q (x,t). The QE states are the eigenstates of the evolution operator U for one period of 
oscillation of the external field T = 2it/u), 

is T 

U(T)^ q {x,t) = exp( —)^ q (x,t), 

which can be defined by (see, for example, Ref.0, p. 385) 

2£ £ t 

i/> q (x, t) = exp( — j-) C*(t)if> n (x) = exp( — ^-)u q (x, t), (7) 

where the functions C%(t) and u q (x, t) are periodic in time, u q (x,t + T) = u q (x,t). 

The coefficients C%(t) are the eigenvectors of the operator U in the representation of the 
Hamiltonian Hq and can be found by diagonalizing the corresponding matrix U n<m . The 
following procedure is one way to obtain the matrix elements^. Let the evolution operator 
U act on the initial state C^°\0) = 5 n ,n ; 



U m , n (T)C^\0) = U m , no (T) = C£°>(T). (8) 

The coefficients C£°>(T) can be computed numerically by integration of Eq. (|4]). They form 
a column in the matrix U m ^ no (T). Repetition of this process for initial states, orthogonal 
to the previous one, C™'(0) = 5 n , n ', n' ^ n , fills the matrix U mjno (T). Diagonalization 
of the matrix U myn (T) yields the eigenvalues e q and the eigenvectors C%. The number of 
the Landau levels N in Eqs. (2), (4) included in our computations is equal to the size of 
the evolution operator matrix U m ^ n {T) and, hence, to the number of its eigenstates — QE 
eigenstates. 

Once the eigenvalues e q and the eigenvectors C% are obtained, then we may write the 
evolution operator for one period U m ^ n (T) in the formlll 

UnAT) = E CM* exp(-te q T/h). (9) 
i 

By raising U n y(T) to degree M and using the orthogonality of the eigenvectors C%, one can 
obtain the evolution operator which propagates the system toward M periods U n>n i (MT) by 

U n , n ,(MT) = E^exp(-z^). (10) 

Given U ntn >(MT), the evolution of any initial state C„(0) can be computed by using 

C n (MT) = E U n>n ,(MT)C n ,(0). (11) 

The expressions ( p!Q[) and ( pTT]) are much more practical for calculations than the integration 
of the set of differential equations (|4]) , especially in the limit t — > oo, because they allow 
us to obtain the state of the system at any time t by a simple summation. There are 
only two dimensionless parameters determining the dynamics of the system, namely the 
dimensionless amplitude of the perturbations Vq = Vq/Huj and the effective Planck constant 
h in the arguments of the matrix elements ([5a]) and (|5b|). This is easy to see if we write Eq.(4) 
in dimensionless form by introducing the dimensionless time r = tuVo- In this form, the 
phases of the oscillating terms are given by — tut/Vq. Thus, the larger the amplitude of the 
wave Vq is, the smaller the frequency of the oscillations becomes and the larger the number 



of the effective terms which participate in the dynamics. Consequently, the parameter Vo 
determines the number of effective terms on the right-hand side of Eq. which can be 
roughly estimated asm > Vo. The same estimation yields the bandwidth of the evolution 
operator matrix @ or (^). If Vo <C 1 the resonant terms with m = ±1 dominate the 
dynamics; they become independent of time since they are being multiplied by sm(ut) or 
cos(ujt). The other terms oscillate fast and can be averaged out. Consideration of only 
these time-independent coefficients constitutes what we call the resonance approximation. 
As was shown in Refs.0El, the quantum resonance approximation in the quasiclassical limit 
in general corresponds to the classical resonance approximation (see e.g. RefH). It is 
necessary to point out that when V > 1, the matrix elements V n>n+m ~ J m {\/2nh) decrease 
quickly with m in the region m > y2nh and this inequality can be treated as an additional 
restriction on the bandwidth of the evolution operator matrix. 

As Vo increases, the number of effective terms increases too. This number was numeri- 
cally determined by requiring relatively small changes in the QE spectrum and in the QE 
eigenfunctions after including an additional term to sum (|j) . The fluctuations did not exceed 
the accuracy of the Runge-Kutta method used for integration of Eq. (Q). The numerically 
determined number of effective terms was found to be of the order of Vo in agreement with 
the speculations presented above. The Runge-Kutta procedure was controlled by the nor- 
malization condition Y^ n \ C n \ 2 — 1; the fluctuations of this value were smaller than 10 -4 . 



III. THE SEPARATRIX QUASIENERGY EIGENSTATES 

In this section we discuss the structure of the QE states which determine the dynamics via 
Eqs. fllPp and (|TT|). We first consider the resonance approximation which will be the starting 
point for the investigation of quantum chaotic effects in next order approximations, ffl The 
equation for the QE eigenstates in the resonance approximation can be obtained by putting 
the QE eigenfunction in the form (7) into Eq. ([|) and keeping only time-independent 
(resonance) terms. Thus, the set of differential-difference equations (^) is transformed into 
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the set of algebraic equations 

EqC^ = V {V n ^ n+ iCl +l + V^n-lC^^), (12) 

where E q = e q /hu is the dimensionless quasienergy. This is an eigenvalue problem for the 
Floquet HamiltonianJll Eq. (12) is similar to the Harper equation in the symmetric gauge 
of the vector potential A with periodic off-diagonal modulation of the matrix elements^ 
In our case the off-diagonal modulation is nonperiodic. The dependence of the matrix 
elements V n ^ n+ \ on the Landau number n is shown in the upper part of Fig. 1. Due to 
oscillations of V n ^ n+ i and Ki,n-i with n, the Floquet Hamiltonian matrix for determining 
the QE eigenstates has a cell structure; the boundaries of the cell are given by the 
zeroes of the Bessel function (c.f. Eq. fl6a])). One can easily show that at small values of E g 
( Vo -C 1) the Floquet Hamiltonian matrix ([12]) can be obtained from the Floquet matrix 
(§). Moreover, as will be shown below, the cell structure is maintained even for very large 
values of Vq (Vo ~ 10); an extremely strong perturbation amplitude is required to destroy 
the cells entirely. The regions where the matrix elements V Utn+ i are small (one such region 
is marked with the rectangle in the plot V n ^ n+ i(n) in Fig. 1 ) can be referred to as "quantum 
separatrices" because in the quasiclassical limit their positions in action I correspond to the 
positions of the classical separatrices in phase space .00 These are given by the zeroes of 
the Bessel function of order 1. 

The structure of the QE eigenf unctions can be understood by characterizing each one 
by its center n q = J2n n \Cn\ 2 an d its dispersion o q = (J2 n { n ~ ^?) 2 |^nl 2 ) ■ The plot of 
n q versus a q in the resonance approximation is shown in the left part of Fig. 2 (a); the 
figures on the right are the Poincare surfaces of sections for the classical system with the 
same parameters. Each point in the plot n q (a q ) corresponds to two QE eigenstates (E q and 
—E q ), due to the symmetry of Eq. ( |T2"D under the transformation 

E q -> -E q , CI - {-l) n Cl (13) 

corresponding to the transformation x — x in Eq. (|7|). It is seen that almost all the QE 
eigenfunctions are divided into groups with practically the same n q and different a q . Each 
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group of states belongs to only one resonant cell because the n q for each QE eigenfunction 
is situated in the center of the cell and n q ± a q does not exceed the size of the corresponding 
cell (the boundaries of the cells in Fig. 2 (a) are marked with arrows). It was shown in 
Ref.0 that the Husimi function^ of a QE state with quasienergy E q > is localized in 
the upper part (x > 0) of the phase space and the Husimi function of a state with — E q is 
situated in the lower part (x < 0). Thus, each row on the left part of Fig. 2 (a) corresponds 
to the two symmetrical classical resonance cells shown in the right part of Fig. 2 (a): the 
first row is associated with two classical cells near the point (x = 0, p = 0); the second row 
with the next symmetrical classical cells; and so on. The number of the QE eigenfunctions 
in an individual cell equals approximately the number of the Landau states in this cell. 

Besides the localized eigenfunctions (arranged in rows in Fig. 2 (a)), there are a small 
number (3-4 percent) of delocalized states which cannot be assigned to any particular 
cell; they are represented by the scattered points. These QE eigenfunctions have large 
dispersions, with a q exceeding the size of one cell. Hereafter we shall characterize the 
"localization length" of the QE functions by a q . It was found numerically that the most 
localized eigenfunctions correspond to the largest QE eigenvalues (in absolute value), while 
the most delocalized states correspond to one of the smallest eigenvalues. Calculations with 
other values of the effective Planck constant h have shown that the number of delocalized QE 
eigenstates increases with increasing h which indicates that, as will be shown below, these 
states are of a pure quantum nature with no classical analogs. A representative delocalized 
eigenfunction, the widest one, marked in Fig. 2 (a) with an arrow, is shown in the lower 
part of Fig. 1, and its Husimi function is plotted in Fig. 3. Note that the eigenfunction as 
well as the Husimi function have their maxima in the regions of the classical separatrices 
(c.f. Fig. 2 (a)). Thus, the delocalized states can be identified as "separatrix eigenstates". 
The Husimi function in Fig. 3 is symmetrical with respect to p — > —p but not symmetrical 
under the transformation x — > —x, which turns the Husimi function corresponding to the 
eigenvalue E q into the one corresponding to — E q (see Eq. (|T3"D). The high peaks in Fig. 
3 near the separatrix line x = result from the slowing of the motion of the classical 
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particle which, in turn, increases the probability of finding the particle in this region. The 
existence of the fully delocalized eigenstates of the matrix U n>m (T) is a very interesting, 
nontrivial feature. Initially ( time M — 0), the QE functions, due to their completeness, 
yield a Kroenecker delta 5 nM > in Eq. fllPp. If n and n! belong to different cells, then only 
a small number of the delocalized QE functions provides the cancellation of terms in Eq. 
(10). The condition of completenes (10) (at M=0) serves as a good check for our numerical 
calculations. At short times (M ~ 1) in the resonance approximation, where Vq <C 1, one 
may take into consideration only the elements along the first off-diagonal of the matrix 
Un tm which are of the order of E q ; the elements in the second off-diagonal will be of the 
order of (E q ) 2 and so on. Let us estimate the width Arii of the maxima of the separatrix 
eigenfunctions, where the index i labels the separatrix number. To this end we approximate 
the matrix element V n ^ n+ i near the separatrices by the linear function (see upper part of Fig. 
1) K. n+i — V^ 0ino+ i + a(n — no). Here n is the Landau state number where the value of 
the matrix element V n ^ n+ i is minimum, and a = dV n0tno+ i/dn. This approximation is valid 
when the separatix region An, is smaller than the total number of the Landau states in the 
i-th cell Hi, Arii -C n,. Under this condition, Eq. ( |E2D takes the form 

E 

TT C n = [Ko,n +l + Ol(n ~ n )]C q n+l + [Ko,n +l + - no - l)]C*_ v (14) 

For the separatrix states, the ratio E q /Vo is of the order of 10 -3 and one can omit the left- 
hand side of Eq. ([14]) ; the minimal matrix element Ki ,n +i near the separatrix is also small 
and may be neglected as well. Under these assumptions, Eq. ( |T4"D gives the following relation 
between the coefficients C%: 

_^r^m-2 1 
°-o+m- m+lm _ 1 --- 2 ^o-l> l A& ; 

Here x is the magnitude of the QE eigenfunction in the maximum and m (here odd) 
is the distance from the maximum (see insert in Fig. 1). The reduced equation ( pT5| ) is 
independent of the parameters and cell number; in consequence, the separatrix width Arii 
should be the same for all cells. This was confirmed by our numerical calculations discussed 
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below. In the quasiclassical limit, when h — > and rii — > oo, the relative width of the 
quantum separatrix An^/n, tends to zero, consistent with the classical dynamics. 

To conclude our discussion on the structure of the QE eigenstates in the resonance 
approximation it is necessary to point out that the numerical results obtained by using two 
different approaches — the Floquet Hamiltonian formalism^ in Eq. flT2|) and the Floquet 
formalism at Vq — > described in Sec. II - - lead to the same results with very high 
accuracy. The latter method is nonperturbative and allows us to obtain the solutions at any 
(not necessarily small) amplitude Vq. In the following discussion this approach will be used 
for investigating quantum chaos in our system. 

In order to incorporate chaos one must increase the perturbation amplitude. The struc- 
ture of the QE eigenfunctions in the presence of weak chaos is shown in Figs. 2 (b), and 
2 (c). Note the qualitative difference from the results obtained for the resonance approx- 
imation (Fig. 2 (a)). Increasing the perturbation amplitude to the value Vq = 6 (Fig. 2 
(b)) gives rise to a change in the separatrix eigenstates: their dispersions a q decrease on 
the average. In contrast, localized eigenfunctions (arranged in rows) are slightly affected by 
the perturbation, resulting in a small splitting of the rows. This is in close agreement with 
the classical behavior; namely, the region around the classical separatrix is first and most 
affected by an increase of Vq. 

Further increase of Vq to the value Vq = 13 (Fig. 2 (c)) leads to an increase in the 
number of delocalized states, and on the average, the localization length grows again. The 
most drastic effect of chaos on the quasienergy eigenfunctions was observed in the region 
of Hilbert space corresponding to completely chaotic motion in classical phase space (the 
first two cells in the right part of Fig. 2 (c)). The chaoticity of the corresponding QE 
eigenfunctions is manifested in the apparent random character of the dependence n q (a q ); 
each QE eigenfunction spans both cells so that they cannot be assigned to any particular one. 
Furthermore, our numerical experiments show that the dependence of on the Landau 
number n in the chaotic area is also very irregular. Note that, as in the classical case, 
perturbation affects the various cells differently: cells with small values of Landau numbers 
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n appear to be more affected than those with larger n. 

In Figs. 2 (b), 2 (c) there are no points corresponding to two eigenf unctions as there 
are in Fig. 2 (a) because of the substantial influence of the nonresonant terms on all the 
QE eigenstates. Recall that in the resonance approximation the system has the symmetry 
defined by Eq. (13). Weak chaos lifts this symmetry and splits the rows on the left-hand side 
of Figs. 2 (b) and 2 (c). The 3rd and 5th rows correspond to mixed phase space dynamics, 
shown on the right hand-side of Figs. 2 (b) and 2 (c). 

IV. QUANTUM DIFFUSION VIA THE SEPARATRICES 

The quantum dynamical manifestations of weak chaos are studied in this section by 
means of the nonperturbative technique based on the Floquet formalism discussed in Sec. II. 
In the regime of weak chaos only a small portion of phase space is chaotic. Correspondingly, 
in the quantum model, only a small number of the eigenstates are affected (see Sec. III). 
However, we expect to detect the influence of weak chaos on the dynamics if we consider 
the diffusion along the separatices. This intuition is based on the fact that the diffusion via 
the separatrices is governed by the separatrix QE eigenstates, which are the ones mainly 
affected by the perturbation. In order to check this, let us consider the evolution of an initial 
state C™°(0) = 5 n>no , which is described by (see Eqs. (0) and (|Tl~D ) 

C7(M) = £ C*C* expH£ 9 M), (16) 
q 

where time is measured in the number of external field periods M. The value of the q-th 
oscillating term is determined by the n-th and no-th amplitudes of the q-th QE eigenfunction 
G\. A transition from no-th to n-th Landau level will occur provided both coefficients C% and 
C%* in Eq. (16) are large. If we consider the transitions between the quantum separatrices, 
then the main contribution to the evolution comes from the separatrix QE eigenfunctions as 
they have their maxima at the various separatrix regions, and the diffusion over the quantum 
separatrices is due to their derealization. The number of separatrix QE eigenfunctions is 
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very small but their effect on the diffusion via the separatrices is crucial because localized 
QE eigenfunctions have minima in the vicinity of the separatrices and do not contribute to 
this process. 

First, let us look at the diffusion over the quantum separatrices in the resonance ap- 
proximation where the system possesses no chaos. To this end, we place an initial state 
C n (0) = S n>no in the separatrix region and follow the evolution of this state for a sufficiently 
long time in order to determine the position in Hilbert space of the quantum particle at time 
t — > oo. The largest characteristic time in the system is t max = 2n/u m i n , where io m in is the 
minimal distance between the effective QE eigenvalues, i.e., the eigenvalues corresponding 
to the eigenstates which constitute the initial state and determine the dynamics. As shown 
above, in our case the effective eigenstates are the separatrix ones. After a time t max the 
dynamics enters into a quasistationary regime and it is convenient to time-average the prob- 
ability in the region t ^> t max in order to eliminate the influence of fluctuations and exclude 
transient effects. 

The result of this procedure is plotted in Figs. 4 (a) and (b) for the same h as in Fig. 2 
and for two different values of Vq; the initial state, marked with a large arrow, was situated 
at the Landau level no = 20 in the second cell near the boundary of the first one. From Fig. 
4 (a) one may note that the time-averaged probability distribution (P n ) in the classically 
inaccessible cells (even in the resonance approximation, where Vq — > 0) is comparable with 
that of the initial (second) cell. It was numerically confirmed that if an initial state C m (0) 
is situated anywhere in the central region of a resonance cell, then only an exponentially 
small part of a wave packet tunnels to the neighboring cells and a logarithmic scale (in the 
probability distribution) is necessary in order to recognize the tunneling phenomenon. El The 
second principal point is the observation that the probability distribution in Fig. 4(a) is 
highest around the boundaries of the quantum resonance cells, being always relatively small 
in the central regions. Thus, anomalously intensive tunneling takes place only between 
the quantum separatrices, and we may refer to this process as "diffusion via the quantum 
separatrices." This effect is reminiscent of the diffusion of a classical particle within the 
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stochastic web. That idea is supported by the Husimi function of the initially 5-like wave 
packet after evolution during a rescaled time r = 3000 (see Fig. 5 (a)). This figure reminds 
us of the web structure in its corresponding classical phase space. However, this effect is 
of a pure quantum nature as the classical particle in the resonance approximation has no 
possibility of penetrating from one cell to another. 

The effect of weak chaos on the diffusion over the separatrices is illustrated in Fig. 4 
(b). Inspection of Figs. 4 (a) and 4 (b) shows that an increase of the amplitude leads to a 
decrease of the diffusion rate. This effect of chaos on quantum diffusion is rather unexpected: 
the increase in the perturbation inhibits tunneling instead of intensifying it. This is a direct 
consequence of the partial localization of the separatrix QE eigenfunctions of Fig. 2 (b) 
caused by the presence of weak chaos in the vicinity of the separatrix. This quantum weak 
chaos effect is manifested in the plot of the Husimi function of the wave packet C n (0) = S njTl0 
after evolution during the time r = 3000 ( see Fig. 5 (b)). Hereafter the time t is measured 
in units of T. The third separatrix is not as clearly defined as in Fig. 5 (a) or in Fig. 
3, being partially destroyed by chaos. The Husimi function looks like the classical density 
distribution within the stochastic web in phase space (see Fig. 5.7 in Ref.i). The 4th and 
5th separatrices are absent and the Husimi function appears to be more localized than in 
the resonance approximation. The structure of the Husimi function within the first two cells 
is irregular, consistent with the more developed (not weak) chaos in this region. 

The phenomenon of localization becomes more evident in Fig. 6 where we show the 
evolution of the squared dispersion a 2 as a function of the rescaled time r = V t. The two 
chosen values of the perturbation amplitude Vo correspond, respectively, to the resonance 
approximation and the regime of weak chaos (the values of h and initial conditions in Fig. 
6 are the same as in Figs. 4 (a), (b)). We see from Fig. 6 that the diffusion over the 
separatrices at Vo = 6 is suppressed in comparison with the resonance case, consistent 
with the dynamical picture of Fig. 4. The data shown in Fig. 6 allow us to estimate 
the effective time of saturation of the probability distribution r max = Vot max . As the wave 
packet in Fig. 6 for Vq = 0.002 spreads over all the separatrices, the effective time in 
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this case is defined by the minimal distance between the separatrix QE eigenvalues; in 
the presence of weak chaos (Vo = 6) T max is mainly determined by the minimal distance 
between the QE eigenvalues within the initial cell, because the probability distribution does 
not evolve to the other cells. The separatrix QE eigenvalues are situated near the center 
of the spectrum where the spectrum is dense and the distances between these levels are 



from the different types of effective eigenstates which determine the dynamics, namely, 



We will now show that the diffusion rate can be characterized by two parameters of the 
probability distribution P n in the neighborhood of each separatrix i. These are the maxima 
of P n , denoted by P iy and its width Arij. Plots like Fig. 4(a) suggest that the distribution 
P n in the neighborhood of the separatrix i, averaged over time, may be approximated by a 
gaussian curve. Then the maxima Pi ( width Arij) is given by the height ( width) of the 
corresponding gaussian. Plots of Pj as a function of the perturbation amplitude V give a 
measure of the diffusion rate. Inspection of Fig. 7 reveals an exponential decrease of the 
diffusion rate up to the value V Q ' ~ 6 (V ' depends on h). This provides further evidence of the 
suppression of quantum diffusion over the separatrices due to weak chaos. Further increase 
of Vq results in a growth of the diffusion rate which can be explained by the average increase 
in the localization length discussed in Sec. III. Such behavior as a function of amplitude 
corresponds to the classical situation. The minimum in the curve V Q ' is the crossover from 
the quantum to the classical diffusion. The further the boundary of the resonant cell is 
from the boundary where the initial state was situated, the smaller the minimum is in the 
corresponding curve because chaos destroys a larger number of the delocalized separatrix 
eigenfunctions which contribute to the diffusion. We speculate that the oscillations of the 
curves in Fig. 7 arise from the fact that the diffusion over the quantum separatrices is 
determined by a small number of (separatrix) QE eigenstates, and turn out to be very 
sensitive to changes in the structure of the QE eigenfunctions. 

Now consider the plots of the widths Arii versus the wave amplitude Vq shown in Figs. 




Thus, a large difference in the values of 



Tmax in the two cases of Fig. 6 arises 



(V = 0.002) > r, 



max 



(V = 6). 
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8 (a), (b), and (c) for the same parameters, initial conditions, and separatrices as in Fig. 7. 
The data, presented in Fig. 8, have the following interesting features: (i) Consistent with 



are approximately the same for all the separatrices and for all values of h (the latter was 
confirmed by our calculations with other values of h). (ii) On the average, the widths An; 
do not change significantly with Vo until the amplitude of the perturbation reaches some 
value Vq, after which the An; begin to grow monotonically. Thus, in the region Vq < Vq, 
classical chaos does not affect the width of the quantum separatrices (but does affect the 
diffusion rate in Fig.. 7). (iii) The threshold V ' in the plot of Anj(Vo) is the same as the 
threshold (position of the minimum) in the plot of Pi(Vo) in Fig. 7. At that point the width 
of the quantum separatrix exceeds the width of the classical stochastic web AHi which can 
be approximated byi 



where rii is the center of the i-th cell and the quantity hn in the quasiclassical limit becomes 
the action I. The results of our calculations using Eq. ( |TTD are presented in Figs. 8 (a), (b), 
and (c) with dashed lines. The discrepancy between the quantum and classical curves may 
presumably be attributed to the approximate character of Eq. fll7|) which is valid only in 
the case of an exponentially thin separatrix. Nonetheless, the trend of the quantum curves 
is qualitatively the same as the classical behavior. 



The numerical data and qualitative analysis presented in this paper allow us to make the 
following conclusions about the nature of quantum weak chaos. In the quantum resonance 
approximation we investigated the resonance structure of Hilbert space and also the new 
phenomenom of quantum diffusion via the separatrix. We remark that this quantum dif- 
fusion has no classical analog because classically, the orbits are confined to resonance cells. 



the above considerations (see Eq. (jI3[)), the widths of the quantum separatrices at Vq < V ' 




(17) 



V. CONCLUSION 
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In quantum mechanics, this diffusion results from tunneling across the quantum separatri- 
ces. Accidental intersection of levels of different cells (associated with the overlap of QE 
functions) with eigenvalues E q <C Vo leads to the formation of delocalized QE eigenfunc- 
tions. In other words, the cell structure of the evolution operator matrix (and the Floquet 
Hamiltonian matrix) gives rise to a long range coupling between states of different cells. 
The dynamical manifestation of this effect is an anomalously large diffusion rate between 
the cells via the quantum separatrices. For sufficiently large values of Vo the nonresonant 
terms may be considered as an effective random perturbation which inhibits the long range 
interaction, thereby localizing the QE eigenstates. 

When the perturbation becomes strong enough diffusion is recovered, but now it is of 
a completely different nature. Namely, for large values of the perturbation we observe an 
analog of the classical diffusion within the stochastic web. This is demonstrated by compar- 
ing the widths of the classical and quantum separatrices. The structure of the quasienergy 
eigenstates which explains the diffusion is also different: in the resonance case they have the 
regular form, which is maintained over several cells of Hilbert space, while at large V such 
states are destroyed and diffusion takes place because the localization lengths of many QE 
eigenstates increase on the average. 

In this paper the evolution operator propagating the system toward an arbitrary number 
of periods of the external field is built in the Hamiltonian H basis. Its eigenstates (QE 
eigenstates) are explored under the condition of resonance in the regime of weak chaos. In 
the new phenomenon of diffusion over the quantum separatrices, it was found that a small 
number of delocalized separatrix QE eigenstates play the dominant role. It was shown that 
weak quantum chaos leads to the localization of the separatrix eigenstates and, hence, to 
suppression of quantum diffusion via the separatrices. At large values of perturbation, V > 
Vq', we have observed a recovery of the diffusion which was associated in the quasiclassical 
limit with the growth of the classical stochastic web in phase space. 

It is necessary to point out that the parameters chosen in our numerical experiments 
correspond to actual experimental situations. Acoustic cyclotron resonance can be observed 
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in a 2D electron gas in semiconductor heterostructures subject to a transverse magnetic 
field, and in the field of a longitudinal sound wave. In order to observe this phenomenon the 
electron relaxation time r p must be large enough to satisfy the inequality uj c t p ^> 1. Under 
this condition, and under the condition of cyclotron resonance, one can choose parameters 
which allow the Fermi level rip to be placed at the boundary (quantum separatix) between 
the first and second cells. That is, the argument of the Bessel function Ji(ka^2np) must 
coincide with the first zero of the Bessel function J\. In order to create this situation in an 
experiment one can choose the following experimental parameters: the sound wave frequency 
should be of the order of 10 GHz, the magnetic field H = 2* 10 3 Oe, the effective electron 
mass m* = 0.7m e , and the electron concentration K = 10 11 cm" 2 . These parameters give: 
h = 0.4, n = n F = 20. The value of the parameter V is determined by the deformation 
div(u)~ ku , where u is the acoustic wave amplitude. Thus, the value V = 10 corresponds 
to the wave deformation ku ~ 10~ 4 10 -5 . When the proposed parameters are realized 
in an experiment we predict that quantum chaos will be manifested as, for example, an 
attenuation of the sound wave. 
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VII. FIGURE CAPTIONS 

FIG. 1. The matrix elements V^n+i in dimensionless units versus the Landau number n (upper 
part); the most delocalized QE eigenfunction in the resonance approximation (lower part). The 
insert amplifies a small portion of the QE eigenfunction, marked with brackets. Here h = 0.37, 
V = 0.002 

FIG. 2. Plot of the dispersions a q versus n q for the QE eigenfunctions with q = 1, 2, . . . , iV 
for h = 0.37, N = 381 and different values of V . (a) V = 0.002, (b) V = 6, (c) V = 13 (left 
hand-side); classical phase space for the same parameters (right hand-side). 

FIG. 3. The Husimi function corresponding to the QE eigenstate plotted in Fig. 1. 

FIG. 4. The time-averaged probability distribution (P n ) versus n, h = 0.37, N = 381, (a) 
Vq = 0.002, (b) Vq = 6. Averaging was performed over 300 times in the region 3 000 < r < 6 000 
(r = Vot; t is measured in units of the external field period T = 2tt/u). The separatrix positions 
are marked with arrows. 

FIG. 5. The Husimi functions of the states for the same parameters and initial conditions as in 
Fig. 4 at (dimensionless) time r = 3000, C n (0) = <5 n! „ , n = 20, (a) V = 0.002, (b) V = 6. 

FIG. 6. The dynamics of the squared dispersion cr 2 (r) for the same parameters and initial 
conditions as in Figs. 4(a) and 4 (b) in units of dimensionless time. 

FIG. 7. The diffusion rate Pi of an initially 5- like wave packet, placed at no = 20, from the 
initial separatrix (first) to the i-th one for i = 3, 4 , 5 as a function of the perturbation amplitude 
Vq. h = 0.37, N = 401. 
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FIG. 8. The widths of the quantum separatrices Arij versus the perturbation amplitude Vq for 
the same parameters, initial conditions and separatrices as in Fig. 7. 
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